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A  distance  propagation  method  is  presented  for  calculating  tortuosity  with  relatively  low  computation 
time  from  three-dimensional  (3D)  tomographic  data.  Moreover,  a  novel  concept  of  tortuosity  distribution 
is  developed  to  provide  a  more  comprehensive  picture  of  inhomogeneous  microstructures  where  tor¬ 
tuosity  depends  on  the  actual  3D  paths.  Instead  of  using  one  single  tortuosity  value,  the  tortuosity 
distribution  both  as  spatial  distribution  map  and  also  statistic  histogram  can  provide  a  more  complete 
description.  The  method,  which  can  be  applied  to  any  porous  medium,  is  tested  against  a  diffusion-based 
tortuosity  calculation  on  two  3D  microstructures:  a  LiCo02  cathode  electrode  of  lithium  ion  battery 
measured  by  x-ray  nano-tomography  and  a  lanthanum  strontium  manganite— yttria-stabilized  zirconia, 
solid  oxide  fuel  cells  cathode  measured  using  focused  ion  beam-scanning  electron  microscopy  serial 
sectioning.  The  present  method  is  shown  to  provide  good-agreement  with  the  effective  diffusion-based 
tortuosity  values. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Tortuosity  is  an  important  parameter  for  characterizing  trans¬ 
port  properties  of  multi-phase  structures,  and  is  of  interest  in  a 
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E-mail  address:  junwang@bnl.gov  (j.  Wang).  terials  [1-4],  catalysts  [5,6]  and  soil  mechanisms  [7-9].  Tortuosity, 
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Fig.  1.  Three  different  voxel  neighboring  definitions:  (A)  city-block,  (B)  quasi- Euclidean  and  (C)  chess-board. 


defined  as  the  ratio  between  the  actual  tortuous  path  length  to  the 
straight-though  path  length,  is  the  parameter  that  embodies  how 
the  specific  phase  structure  influences  transport.  As  three  dimen¬ 
sional  (3D)  structural  measurements  have  become  increasingly 
available,  it  has  become  useful  to  calculate  tortuosity  values  directly 
from  measured  3D  structures.  One  example  is  for  characterizing, 
modeling,  and  designing  energy  storage  and  conversion  materials 
such  as  batteries  and  fuel  cells,  where  transport  within  the  elec¬ 
trodes  is  directly  relevant  to  electrochemical  properties. 

The  importance  of  tortuosity  for  battery  electrodes  has  been 
discussed  previously  [1,2,10-17  .  Given  that  active  electrode  ma¬ 
terials  in  Li-ion  batteries  make  up  typically  50-70%  of  electrode 
volume,  and  that  the  remaining  volume  is  shared  between  elec¬ 
trolyte,  binder,  and  carbon  phases,  Li-ion  transport  in  the  electro¬ 
lyte  phase  may  limit  battery  capacity,  compared  with  the  electronic 


conductivity  and  charge-transfer  [18,19],  in  some  cases.  Attempts  to 
increase  the  density  of  active  materials  or  increase  electrode 
thickness,  in  order  to  increase  the  energy  density,  will  ultimately  be 
limited  by  maintaining  the  ion  transport  [1].  Electrolyte-phase 
tortuosity  is  thus  a  key  factor  in  design  of  electrodes.  Advanced 
processing  approaches  have  been  used  to  design  battery  electrode 
morphology  and  architecture  to  improve  power  performance  via 
minimizing  the  electrode  tortuosity  and  therefore  to  maximize 
power  with  a  given  porosity  and  thickness  of  the  electrodes  [1  .  In 
addition,  others  have  used  tortuosity  as  an  effective  geometric 
parameter  in  models  to  describe  ionic  transport  in  the  porous 
electrode  [2,12,20,21].  Similarly,  the  importance  of  tortuosity  has 
also  been  discussed  in  many  literature  regarding  energy  conversion 
materials  such  as  proton  exchange  membrane  fuel  cells  and  solid 
oxide  fuel  cells  (SOFC)  [3,4,22-27].  In  SOFC  electrodes,  tortuosity  of 
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Fig.  2.  Four  different  methods  to  extract  the  tortuosity  from  the  distance  curve.  (A)  Scale:  the  curve  is  fitted  with:  y  =  mx  +  b  where  x  is  the  straight  distance,  y  is  the  path  length 
and  the  tortuosity  is  m.  (B)  End:  tortuosity  is  defined  by  average  path  length  at  the  end  plane  divided  by  the  straight  distance.  (C)  Average:  tortuosity  is  defined  by  the  average 
tortuosity  from  each  plane;  tortuosity  =  1/n  Ya=\  (D)  Slope:  tortuosity  is  extracted  by  average  of  all  local  slopes.  The  local  slopes  are  obtained  from  first  fitting  the  curve  first 
with  a  polynomial  with  degree  of  five  and  then  calculating  the  derivative  at  each  x  (straight  distance)  point. 
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Fig.  3.  Cropped  reconstructed  3D  structure  of  the  (A)  LiCo02  sample  imaged  by  x-ray  nano-tomography.  (B)  70/30  LSM-YSZ  sample  and  (C)  50/50  LSM-YSZ  samples  imaged  by 
FIB-SEM  serial  sectioning. 


the  pore  and  ionically-conducting  phases  may  both  be  important. 
Vapor  transport  in  the  pore  phase  is  directly  related  to  concentra¬ 
tion  polarization  [28],  while  oxygen  ion  transport  in  ionic- 
conductor  phases  can  limit  the  effective  active  thickness  in  the 
electrode  [29]. 

Two  different  types  of  tortuosity  calculation  have  been  reported 
previously.  One  is  the  effective  diffusivity  method,  based  on 
simulating  the  diffusive  or  conductive  transport  across  one  phase  in 
a  measured  3D  microstructure;  dividing  by  the  transport  flux  for 
straight-through  paths  yields  the  tortuosity  factor  [10,23,30].  The 
other  method  is  a  geometrical  calculation  where  the  lengths  of  the 
tortuous  paths  in  the  3D  microstructures  are  measured;  the 
resulting  path  length  divided  by  the  straight-through  path  length  is 
the  tortuosity  3,31].  Path  length  calculation  has  been  done  in 
various  ways,  including  distance  prorogation  [32]  and  shortest  path 
search  3].  While  the  geometrical  method  is  computationally  much 
less  intensive  compared  with  the  effective  diffusivity  method,  it  is 
not  clear  how  well  the  tortuosity  calculated  from  geometrical 
methods  represents  the  diffusion  behavior.  No  direct  comparison 
has  been  reported  yet  to  address  the  difference  or  similarity  be¬ 
tween  these  methods.  Furthermore,  the  prior  calculations  yield 
only  one  tortuosity  value  to  represent  inhomogeneous  structures, 
and  hence  do  not  provide  information  on  the  tortuosity  variation 
within  the  structure  [3,33-36]. 

Here  we  describe  a  path-length  method  to  characterize  the 
tortuosity  of  phases  in  various  complex  3D  microstructures,  and 
compare  the  resulting  values  with  effective  diffusivity  method 
calculations  on  the  same  microstructures.  We  further  developed  a 
concept  of  tortuosity  distribution  by  visualizing  a  spatial  distribu¬ 
tion  of  the  tortuosity,  and  also  quantifying  the  tortuosity  variation 
with  histograms  to  overcome  this  major  drawback.  Various  mi¬ 
crostructures  with  a  range  of  structural  characteristics  and  tortu¬ 
osity  values  are  considered. 


2.  Method 

2.1.  Background 

Tortuosity  is  defined  as  the  ratio  between  the  actual  path  length 
to  the  straight  distance.  For  a  porous  medium  with  Leff  as  the 
effective  actual  path  length  from  point  A  to  point  B,  and  the  L  as  the 
straight  distance  (Euclidean  distance),  the  strict  geometrical  defi¬ 
nition  of  its  tortuosity  (t)  is: 


Geometrical  method  directly  uses  this  definition  to  calculate 
tortuosity  by  measuring  Leff  and  I. 

Historically,  this  geometry  factor  was  introduced  by  solving  the 
one-dimensional  diffusion  problem,  Fields  first  law  of  diffusion: 


where  J  is  the  flux  of  diffusing  species,  D  is  its  diffusion  coefficient,  c 
is  its  concentration  and  x  is  the  position  along  the  diffusion  path.  To 
relate  diffusion  in  porous  media  to  equation  (2),  it  was  assumed 
that  the  transport  media  are  cylindrical  capillaries  [37].  This  leads 
to  the  relationship  between  the  diffusion  coefficient  D  in  a  non- 
tortuous  path  and  the  effective  diffusion  coefficient  Deff  in  a 
tortuous  path: 


Deff 


(3) 
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Fig.  4.  Three-dimensional  (3D)  distance  map  based  on  geometric  prorogation.  The  neighboring  definition  used  here  in  the  figure  is  city-block.  (A)  propagating  along  positive  x 
direction,  (B)  propagating  along  negative  x  direction,  (C)  propagating  along  positive  y  direction  and  (D)  propagating  along  positive  z  direction.  The  3D  distance  maps  for  the  negative 
y  and  z  directions  are  not  shown  here  for  brevity. 


where  £  is  the  volume  fraction  of  the  pores.  The  geometry  factor 
tortuosity  (t)  was  then  defined  as  Leff/L  shown  in  equation  (1).  The 
square  of  tortuosity  is  therefore  called  tortuosity  factor  (T),  defined 
as: 


Though  it  is  common  to  see  the  mix  use  of  the  tortuosity  (t)  and 
the  tortuosity  factor  (T)  in  the  literature,  extra  caution  should  be 
taken  to  avoid  this  confusion.  Tortuosity  calculation  based  on 
effective  diffusivity  method  can  then  be  performed  using  equation 
(3)  or  its  equivalent. 

2.2.  Approach 

In  this  paper,  both  Li-ion  battery  and  SOFC  electrode 
structures  were  used  to  test  the  accuracy  of  the  present  calculation. 
First,  a  porous  LiCo02  structure,  cathode  electrode  of  Li-ion  bat¬ 
tery,  [38]  was  imaged  by  x-ray  nano-tomography  [39].  Second, 
two  lanthanum  strontium  manganite-yttria-stabilized  zirconia 
(LSM-YSZ),  SOFC  cathodes,  with  LSM:YSZ  weight  ratios  of  =  70:30 


(70/30)  and  50:50  (50/50),  were  imaged  and  reconstructed  using 
focused  ion  beam-scanning  electron  microscopy  (FIB-SEM)  serial 
sectioning  [29].  Both  the  geometrical  method  and  effective  diffusion 
methods  were  applied  to  calculate  the  tortuosities  of  these  3D 
structures.  The  following  phases  were  in  the  interest:  the  electrolyte 
(pore)  phase  in  LiCoC^,  and  the  pore  and  YSZ  phases  in  the  LSM-YSZ. 

The  detailed  description  for  the  effective  diffusion  method  can 
be  found  elsewhere  [22].  The  geometrical  method  used  here  to 
calculate  tortuosity  was  based  on  directional  distance  map  propa¬ 
gation.  As  it  is  important  to  calculate  tortuosity  along  different  di¬ 
rections  in  the  microstructure  to  reveal  the  anisotropy  [32],  the 
tortuosity  was  calculated  in  three  orthogonal  directions  (x,  y  and  z). 
Within  the  phase  of  interest,  only  the  interconnected  portions  that 
can  contribute  to  transport  were  included  in  the  calculation. 

Flere  we  elaborate  the  details  of  the  algorithm  by  describing  the 
steps  to  calculate  the  tortuosity  along  the  x  direction.  The  tortuosity 
along  two  other  orthogonal  directions  y  and  z  were  calculated  using 
the  same  procedure.  First,  all  of  the  phase  of  interest  located  in  the 
starting  plane,  ‘seed  plane’,  are  labeled  with  distance  of  unity.  Sec¬ 
ond,  the  neighboring  voxels  of  the  phase  of  interest  are  labeled  with 
their  effective  distance,  and  this  procedure  is  repeated  directionally 
through  the  phase  network  until  all  the  voxels  is  labeled  with  a 
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Fig.  5.  Distance  curves  of  LiCo02  show  the  actual  path  length  vs.  straight-through  distance  for  all  three  orthogonal  directions,  prorogating  from  (A)  positive  and  (B)  negative 
directions. 


distance.  Depending  on  the  neighboring  voxel  definition,  the 
effective  distance  may  vary.  Here  we  adapted  three  different 
neighboring  definitions:  1)  city  block,  2)  quasi-Euclidean  and  3) 
chess  board.  A  schematic  is  shown  in  Fig.  1  to  illustrate  the  concept  of 
different  neighboring  definitions  and  how  it  affects  the  distance 
propagation.  Third,  the  so-called  ‘distance  curve’,  i.e.,  a  plot  of  the 
effective  lengths  versus  the  straight-through  distance  was  obtained. 

Finally,  the  tortuosity  was  calculated  from  the  distance  curve. 
We  implemented  four  different  methods  to  extract  the  tortuosity. 
The  first  method  is  by  fitting  the  curve  with  one  single  slope  as  used 
previously  in  the  literature  [31  ].  Two  problems  arise  from  defining  a 
single  fitted  slope  as  tortuosity.  First,  there  is  a  possibility  to  obtain 
a  slope  of  less  than  unity,  which  is  physically  impossible  as  an 
actual  path  cannot  be  less  than  the  straight  distance.  Secondly,  it 
only  represents  an  overall  trend  and  the  step-determining  slow 
diffusion  structures  are  not  reflected  in  the  result.  To  resolve  these 
problems,  three  more  methods  are  proposed  and  tested  here.  The 
illustration  of  the  different  methods  to  extract  tortuosity  from  the 
distance  curve  is  shown  in  Fig.  2.  Other  methods  are  named  with 
End,  Average  and  Slope. 

Although  our  method  is  based  upon  the  distance  propagation, 
we  have  modified  it  in  the  following  aspects.  First,  the  distance 
propagation  was  performed  with  the  entire  volume  instead  of  from 
the  center.  As  previously  used  in  the  literature,  the  seed  plane  is 
defined  as  the  center  plane  of  the  sample  and  therefore  the 


algorithm  effectively  measures  two  sub-volumes,  which  is  different 
from  measuring  from  the  entire  volume.  Here,  the  distance  prop¬ 
agation  is  performed  separately  along  the  positive  and  negative 
directions  and  then  the  average  is  also  calculated.  Secondly,  we 
included  quasi-Euclidean  distance  to  improve  the  accuracy  in 
measuring  the  distance  between  voxels.  Previously  by  using  city- 
block  definition,  the  tortuosity  tends  to  be  over-estimated  and  by 
using  chess-board  definition,  the  tortuosity  tends  to  be  under¬ 
estimated.  Last,  different  methods  were  used  to  extract  the  tortu¬ 
osity  instead  of  only  using  a  single  slope  as  its  limitation  pointed 
out  previously. 

With  the  distance  map  obtained  via  distance  propagation,  the 
following  distributions  can  also  be  obtained  directly:  1 )  histogram 
of  tortuosity  at  a  given  plane  and  2)  spatial  distribution  of  the 
tortuosity,  viewed  as  a  tortuosity  map.  This  additional  information 
provides  a  comprehensive  picture  to  characterize  the  tortuosity 
and  to  facilitate  the  understanding  of  the  transport  phenomena  in 
the  porous  media. 

3.  Results  and  discussion 

3.1.  Propagation  of  distance  map 

The  reconstructed  3D  microstructures  of  the  LiCo02  and  LSM- 
YSZ  electrodes  are  shown  in  Fig.  3.  Fig.  4  shows  the  3D  distance 


Fig.  6.  Comparison  of  the  geometry  method  vs.  diffusion  method  for  LiCo02  pore  phase.  The  geometry  method  here  used  the  city-block  as  the  neighboring  definition. 
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Fig.  7.  Comparison  of  the  geometry  method  vs.  diffusion  method  for  LSM-YSZ  pore 
and  YSZ  phases,  both  70/30  and  50/50  samples.  The  geometry  method  here  using  city- 
block  as  the  neighboring  definition  and  average  as  the  calculation  method. 

maps  based  on  geometric  distance  prorogation  along  the  pores  of 
LiCo02  for  all  three  orthogonal  directions.  For  each  direction,  the 
path  distance  maps  along  both  positive  and  negative  directions 
were  calculated.  For  instance,  Fig.  4(A,  B)  shows  the  distance 
propagation  along  the  positive  and  negative  direction  of  the  x  axis 
respectively.  The  similar  distance  maps  were  also  obtained  for  the  y 
and  z  directions  but  only  the  distance  maps  along  the  positive  di¬ 
rections  are  shown  here  for  brevity.  The  distance  maps  were  also 
obtained  for  the  LSM-YSZ  sample  and  are  not  shown  here  for 
brevity. 

After  establishing  the  distance  map,  the  average  actual  path 
length  at  each  plane  for  a  fixed  straight  distance  along  x,  y  and  z  can 
then  be  obtained.  The  results  using  the  city-block  neighboring 
definition  are  shown  in  Fig.  5.  Similarly,  the  distance  curve  for  the 
other  two  voxel  neighboring  definitions  were  also  calculated  but 
are  not  shown  here  for  brevity. 

3.2.  Comparison  between  the  geometry  method  and  the  diffusion 
method 

The  comparison  of  the  geometrical  and  effective  diffusion 
method  results  for  LiCo02  is  shown  in  Fig.  6.  The  geometry  method 
shown  here  is  using  the  city-block  method  as  the  neighboring 
definition,  which  is  found  to  be  closest  to  the  tortuosity  values 
calculated  from  the  diffusion  method.  The  effects  of  the 


neighboring  definitions  are  discussed  in  Section  3.3.  Among  the 
four  different  definitions  to  extract  the  tortuosity  in  the  geometry 
method,  the  definition  that  produced  the  results  closest  to  the 
diffusion  method  was  the  average  method.  This  indicates  that  the 
diffusion  method  has  a  similar  nature  as  the  geometry  method 
when  taking  into  account  the  effect  of  each  of  the  individual  planes 
and  defining  the  total  tortuosity  as  the  average  of  the  tortuosity 
from  each  plane. 

The  tortuosity  results  of  the  LSM-YSZ  samples  from  both  the 
geometry  and  diffusion  methods  are  shown  in  Fig.  7.  These  two 
methods  show  good  agreement  overall  with  the  trend  when 
comparing  the  tortuosity  between  different  phases  or  between 
different  directions.  For  the  pore  phase  where  the  tortuosity  is 
lower,  these  two  methods  show  better  agreement  in  their  values. 

3.3.  Effects  of  different  parameters  in  the  geometry  method 

As  shown  previously,  the  four  different  methods  to  calculate 
tortuosity  from  the  distance  curve  lead  to  slightly  different  tortu¬ 
osity  results  since  they  show  different  aspects  when  representing 
tortuosity.  In  addition  to  this  parameter,  other  parameters  that 
might  affect  the  calculations  include:  distance  propagation  direc¬ 
tion,  sampling  volume  size,  and  neighboring  voxel  definition, 
which  are  discussed  in  the  following. 

For  a  structure  that  is  fully  connected  in  3D,  meaning  that  the 
distance  propagation  can  be  performed  completely  from  both 
positive  and  negative  directions,  the  difference  that  results  from 
the  direction  definition  is  generally  less  than  2%,  independent  of 
other  factors  such  as  the  neighboring  voxel  definition  and  also  of 
the  sampling  volume  size. 

The  sampling  volume  size  affects  the  tortuosity  simply  due  to 
the  inhomogeneity  in  the  studied  media.  Figure  8(A)  shows  the 
relationship  between  the  tortuosity  of  the  pores  versus  the  sample 
volume  size  for  the  LiCoCb  sample.  The  variation  of  the  tortuosity  as 
a  function  of  the  sampling  volume  is  therefore  an  indication  of  the 
inhomogeneity  of  the  sample. 

The  neighboring  voxel  definition  affects  the  tortuosity  in  the 
calculation  of  the  actual  path  length.  A  comparison  between  the 
three  different  neighboring  definitions  is  shown  in  Fig.  8(B).  By 
varying  the  neighboring  voxel  definition,  the  physical  connectivity 
between  two  nearby  voxels  varies.  It  is  natural  that  the  most  strict 
neighboring  definition  city-block  results  in  the  largest  tortuosity 
since  only  the  face-adjunct  voxels  are  considered  neighbors.  It 
follows  by  the  results  from  quasi-Euclidean  which  is  a  more  neutral 
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Fig.  8.  (A)  Tortuosity  vs.  sampling  volume  size  for  LiCo02  sample.  Here  the  sampling  volume  is  the  length  of  the  cubic  volume  that  is  being  sampled.  Here  the  sampling  volume  is 
the  length  of  the  cubic  volume  that  is  being  sampled.  For  instance,  the  label  100  refers  to  a  volume  with  100  x  100  x  100  voxels.  (B)  The  tortuosity  results  depend  on  the  voxel 
neighboring  definition  -  city-block,  quasi-Euclidean  and  chess-board. 
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Fig.  9.  “Tortuosity  distribution”  of  the  LiCo02  sample  shown  as  (A-C)  spatial  distribution  maps  and  also  (D-F)  histograms  for  all  three  orthogonal  directions  along  the  positive 
direction. 


definition  of  the  neighboring  voxel  distance.  The  chess-board 
definition  is  found  to  tend  to  underestimate  the  tortuosity  signifi¬ 
cantly  and  should  not  be  used. 

3.4.  Tortuosity  distribution  characterized  by  spatial  distribution 
map  and  histogram 

In  order  to  overcome  the  limitation  of  using  one  single  tortu¬ 
osity  value  to  represent  the  whole  structure,  the  “tortuosity  dis¬ 
tribution”  both  as  a  spatial  distribution  map  and  also  a  histogram 
can  be  obtained  from  the  distance  propagation  map  shown  in  Fig.  4. 
The  tortuosity  distribution  shown  as  spatial  distribution  maps  for 
all  three  orthogonal  directions  of  UC0O2  pore  phase  are  presented 
in  Fig.  9(A— C).  The  significant  variation  can  be  observed  clearly 
with  some  regions  of  the  structure  has  tortuosity  as  high  as  almost 
2.5  and  other  regions  of  the  structure  as  low  as  almost  unity.  This 
spatial  variation  was  lost  when  averaging  the  distance  along  a 
specific  distance  plane  and  cannot  be  characterized  by  using  one 
single  tortuosity  value.  This  tortuosity  variation  can  also  be  char¬ 
acterized  by  histograms  as  shown  in  Fig.  9(D-F).  A  single  value 
tortuosity  only  shows  the  mean  value  of  the  tortuosity  distribution. 

4.  Conclusions 

A  geometrical  tortuosity  calculation  method  was  tested  on  three 
measured  3D  structures:  a  UC0O2  battery  electrode  and  two  3D 
LSM-YSZ  structures.  The  comparison  between  the  geometrical  and 
diffusivity  methods  was  discussed,  and  the  specific  geometrical 
method  yielding  the  best  agreement  was  identified.  The  method 
and  the  concept  of  tortuosity  distribution  based  on  large  volume  of 
3D  structure  using  x-ray  nano-tomography  opens  up  a  way  of  ac¬ 
curate  and  efficient  characterization  of  the  tortuosity.  Although  the 
present  examples  are  energy  materials,  the  method  is  broadly 
applicable  to  any  porous  or  multi-phase  media  studies.  The  effec¬ 
tive  diffusivity  can  be  derived  for  a  given  species  in  a  given  media, 
once  the  tortuosity  is  known.  The  tortuosity  can  also  be  used  as  a 
critical  parameter  in  modeling  and  designing  of  energy  storage  and 
conversion  materials. 
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